Reflection of two-gap nature in penetration depth measurements of MgB2 film 
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The magnetic penetration depth, A(r), in the basal plane of a magnesium diboride (MgB2) film 
was measured using a two-coil mutual inductance technique at 50 kHz. This film has Tc — 38 K, 
ATc < 1 K, and A(0) ~ 1500 A. At low temperatures, A~'^(r) shows a clear exponential temperature 
dependence, indicating s-wave superconducting order parameter symmetry. However, the data 
are not quantitatively well described by theory assuming a single gap. From the data fit by the 
full BCS calculation assuming a double gap, the values of the two distinct gaps were obtained: 
As(0) = 2.61 ± 0.41 meV and Al(0) = 6.50 ± 0.33 meV. The contributions of the smaU and the 
large gaps to the total superfiuid density at T = were estimated to be 21% and 79%, respectively. 
Finally, we consider the effect of gap anisotropy on the penetration depth measurements, and find 
that the gap anisotropy does not play a significant role in determining the temperature dependence 
of the penetration depth. 



I. INTRODUCTION 



Since the discovery of superconductivity in metallic MgB2, abundant research has been carried out to elucidate 
its basic mechanism. MgB2 has some notable features which contrast with cuprate superconductors. First, electron 
coupling is mediated by phonons, as indicated by the observation of a prominent isotope effect. Q However, the 
transition temperature, Tc ~ 39K, might be somewhat higher than the theoretical prediction assuming a conventional 
phonon mechanism. Recently, it was suggested that anisotropy in electron-phonon couplings plays a significant 
role in the unusually high transition temperature. iQ,^ Secondly, many theoretical jJ-Q] and experimental works 
[p|-pO[ suggest that MgB2 has two separate gaps and that the symmetry of each gap is s-wave with a substantial gap 
anisotropy. The Fermi surface of MgB2 consists of two nearly cylindrical (2D) sheets and two tubular networks (3D). 
PH While the value of the gap associated with 3D sheets is in the range 0.5 < A{0) /ksTc < 0.95, for 2D networks 
the value is about 1.8 < A{0) /ksTc < 2.2. This double-gap structure and its anisotropic nature would be expected 
to play an important role in the physical properties of the compound. 

In the early stages of MgB2 study, a number of groups claimed unconventional superconductivity or s-wave 
order parameter symmetry p3|-p5| from penetration depth measurements on various forms of samples. For example, 
Prozorov et al. |25| showed a clear exponential behavior of A(T) measured using a microwave technique for MgB2 
wires. On the other hand, Pronin et al. p^ ] reported the penetration depth measurements of c-axis oriented films 
with Tc ~ 32 K. They claimed that the penetration depth shows temperature dependence rather than exponential 
behavior at low temperatures, and suggested a strong gap anisotropy or the existence of nodes in the gap as its origin. 
At present, some of the results of penetration depth measurements and their interpretations are still controversial. 

In this work, we measure the magnetic penetration depth of a high-quality MgB2 film via a mutual-inductance 
technique. At low temperatures, the superfiuid density, ps oc A~^(T), shows a clear exponential temperature de- 
pendence. However, the model assuming a single gap does not describe the experimental results. A full calculation 
within the BCS framework assuming existence of a double gap, successfully describes our data. The values of the two 
gaps obtained from fits to the data are consistent with previous reports. Also, the contributions of each gap to the 
superfiuid density are deduced from the analysis. 

Initially, we compare the data with the low-temperature expansion formula for A^^(T) assuming one or two gaps 
to get a rough estimate of the gap values. Then, we refine the values from a full calculation of the penetration depth. 
Finally, we consider how the penetration depth is affected by gap anisotropy. 
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II. EXPERIMENTAL ASPECTS 



The MgB2 thin fihii was fabricated using a two-step method; the detailed process is described elsewhere. First, 
an amorphous boron thin film was deposited on a (1 1 2) AI2O3 substrate of 1 x 1 cm^ at room temperature by 
pulsed laser. Then, the boron thin film was put into a Nb tube with high purity Mg metal (99.9%), and the Nb tube 
was then sealed using an arc furnace in an Ar atmosphere. Finally, the heat treatment was carried out at 900°C for 
10 to 30 minutes in an evacuated quartz ampoule, which was sealed under high vacuum. The film thickness is 0.3 
/im, confirmed by scanning electron microscopy. X-ray diffraction patterns indicated that the MgB2 thin film has a 
highly c-axis-oriented crystal structure normal to the substrate surface; no impurity phase is observed. 

The penetration depth, A(r), was measured using a two-coil mutual inductance technique described in detail 
elsewhere. |^,^ The MgB2 film is centered between drive and pick-up coils with diameter of 1 mm. The inset of 
Fig. |l| schematically illustrates the measurement configuration. A current (1 mA< Id < 30 mA) at 50 kHz in the drive 
coil induces screening currents in the film. The net magnetic field from the drive coil and the induced current in film 
are measured as a voltage across the pick-up coil. Because the coils are much smaller than the film, the applied field 
is concentrated near the center of the film and demagnetizing effects at the film perimeter are not relevant. All data 
presented here are taken in the linear response regime. Figure |l| shows representative mutual inductance, M{T), data 
measured with Id — 30 mA for T < 25 K and Id — 2 mA for T > 25 K. The mutual inductance technique enables us 
to extract absolute values as well as temperature dependence of the penetration depth from the mutual inductance 
data. 

The procedure to extract A~^(T) from M{T) is the following: First, a constant background (zero position) due to 
stray couplings between coils is subtracted from raw data. This constant background can be estimated from measuring 
at r = 4.2 K the mutual inductance of Pb foil with identical shape and area as substrate using the same measurement 
probe. In this background measurement, the magnetic penetration depth of Pb is so small compared to foil thickness 
that no magnetic field goes through film. After the subtraction of the background, the data is normalized to the value 
of mutual inductance at T = 50 K (initial position) . The normalization removes uncertainties associated with amplifier 
gains and nonideal aspects of the coil windings. The subtracted and normalized mutual inductance is converted to 
complex conductivity, a = ai — ia^ , where cti and cr2 are real and imaginary parts of the conductivity. Finally, the 
penetration depth is determined from the imaginary part of conductivity via the relationship a2 = 1/ ^a^X"^, where 
/iQ is the magnetic permeability of vacuum and oj is the frequency of drive current. The accuracy of is limited by 
10% uncertainty in film thickness. However, the temperature dependence of is unaffected by the uncertainty. The 
other inset of Fig. |l| displays A~^(T) curves deduced from M(T) measured at two different Id levels. Although the 
signal-to- noise ratio of upper curve (/^ = 10 mA for T < 25 K) is smaller than that of lower one {Id = 30 mA for T < 
25 K), two curves do not show any quantitative difference. (Upper curve was shifted by 10 /j,m~^ for comparison.) 



III. RESULTS AND DISCUSSION 



A. Theoretical description of data assuming a single gap 

Figure | shows X-'^{T) at temperatures below 15 K. The value of X^'^iT) at T ~ 1.3 K is about 43 /im , which 
corresponds to A c± 150 nm. To examine the temperature dependence of A^^(T) at low temperatures, we fit first ~ 5% 
drop in A^^(r) to an exponential-type function A^^(T) ~ 1 — cexp(— Z3/T), where c are D are adjustable parameters. 
As presented by thin solid line in the figure, the fit is reasonably good. On the other hand, when we fit the data in 
the same temperature region to a quadratic form A~^(r) 1 — (T/Tq)^ as in Ref. [^,|6|, the fit deviates significantly 
from the data (dashed line in the same figure). The values for exponential and quadratic fits are 3.96 x 10~^ and 
5.13 X 10~*, respectively. 

The exponential T-dependence of X~^{T) at low temperatures can be regarded to refiect s-wave order-parameter 
symmetry in this compound. For clean BCS-type superconductors A^^(T) is given by 

(1) 



A-2(G) V dE 

where A is the superconducting energy gap, f{E) = [exp(— iJ/fcsT") + 1]^^ is the Fermi distribution function, and 
D{E) = E/{E^ — A^)^/^, is the quasiparticle density of states (DOS). This equation can be expanded at low temper- 
atures, where A is nearly constant, in the following way |]3l[|: 



A-'(r) / 2^A(o) y/' 

A-2(0) - [ keT J 



exp(-A(0)/fcBT), (2) 
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where A(0) is the energy gap at zero temperature. Thin soUd hne in Fig. ^ represents Eq.(||) fitted to data. The 
comparison between data and theory is restricted to the first ~ 5% drop in A^^(T), where the T-dependence of A 
is not significant. This comparison yields A(0) = 4.29 meV [A(0)/fcBT'c = 1-31] and \^^{0) = 43.2 /im^^ xhe gap 
value is significantly smaller than the BCS weak coupling limit A{Q) /kBTc ~ 1.76. Using the gap value deduced 
above, A^^(r) in the whole temperature region below Tc can be obtained by a full BCS calculation using Eq. (|^). 
If the above one-gap fit is valid, the full calculation is expected to describe the experimental A~^(T) for the entire 
temperature region below Tc- Thin solid line in the inset of Fig. ^ represents this full calculation. The curve does 
not give a correct description of the data at high temperatures. While the data show negative curvature at high 
temperatures, the theoretical line shows weakly positive curvature. 



B. Theoretical description of data assuming a double gap 



A number of experimental and theoretical groups have proposed the existence of two gaps in the DOS of MgB2. 
The larger gap belongs to the quasi-2D Fermi surface derived from B-B (tr) bonds, and the smaller gap belongs to the 
quasi-3D Fermi surface derived from B-Mg-B (tt) bonds. We model this two-gap nature by writing D{E) as a sum of 
two BCS DOS. Thus, low-temperature expansion of A^^(r) can be expressed by 



A-2(T) 



1 - 



2^As(0)V^' 



+AZ'(0) 



27rAL(0) 
knT 



exp(-As(0)/fcsT) 

1/2 

exp(-Ai(0)/fcsr) 



where the A5(0) [Al(0)] and A|^^(0) [A£^(0)] are the value of the small (large) gap and the contribution of small 
(large) gap to total superfluid density (oc A^^), respectively. For comparison with data, it is more convenient to 
convert the above equation to the following form 



A~^(T) 
A-2(0) 



1-ci 



27rAs(0) 



-C2 



keT 

/^2 7rAL(0) 
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exp(-AL(0)/fcBT), 



(3) 



where X-^{0) = XfiO) + A;^^(0), ci = Xf{0)/X-^{0), and C2 = (1 - ci). 

At very low temperatures, the change of superfluid density with temperature, i.e., the quasiparticle excitation, is 
dominated by the small gap. In other words, the role of the large gap in A^^(r) is relevant at higher temperatures. 
Thus, we extend the fitting region up to ~ 15% drop in A~^(T). 15% drop in A^^(T) corresponds to about 10% drop 
in A(T)/A(0) in the BCS weak coupling limit. Accordingly, in this fit about 10% error due to change of A can be 
expected. Figure § shows the comparison of Eq. (§) with data. From this fit, we obtains two distinct gap values, 
As(0) = 2.57 meV and Al(0) = 5.82 meV, corresponding to As{0)/kBTc = 0.79 and AL(0)/fcsTc = 1.78. In the 
case of the small gap, the value is consistent with previous reports. But the large gap value is somewhat smaller than 
those in the literatures (Table |). 

The success of the fit motivates a full calculation in the extended temperature range for a more precise description 
of the data. We assume isotropic s-wave gaps on the two pieces of Fermi surface, and perform a full calculation of 
A~^(T) according to 



A-^(r) 

A-2(0) 



= 1-2 



ci 



(4) 



where Ci is adjustable parameter which determines the contribution of the small gap to the superfluid density and 

C2 = (1 - Ci). 

Figure ^ shows our attempt to fit the data using Eq.(Q). Except near Tc, the theoretical line gives a good fit to 
the data. From this, we obtain the gap values A5(0) = 2.61 ± 0.41 meV and Ai(0) = 6.50 ± 0.33 meV. These are 
fairly consistent with previous reports (Table I). Also, the contributions of each gap to A^^(O), i.e., C2 = 0.79 ± 0.06 
is deduced. The inset of Fig. || shows theoretical A^^(T) curves, where the contributions of each gap are separately 
plotted. 
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In the above analysis, we described the A~^(T) theoretically assuming the gaps (A^ and A5) being isotropic on the 
Fermi surfaces. According to a recent theoretical calculation, the values of the small and large gaps are distributed 
in the range of 1 meV < A5 < 3 meV and 6.5 meV < A^ < 7.5 on the Fermi surfaces. 1^ Here we suppose two 
phenomenological models for gap distribution. In the first model, the gap is distributed uniformly around the average 
value of gap, Aq. In the alternative model, we assume the normal (Gaussian) distribution of the gap. Using these 
models, we calculate theoretical curves of A^^(T) assuming ~ ±25% variation of gap around Aq on the Fermi surface. 
The gap distribution of ±25% is sufficient to account for real gap anisotropy in MgB2. The calculations reveal 
that the change in A~^(T) due to the gap distribution is not significant. In fact, the maximum change in due 
to the gap anisotropy is only about 2% in the case of the uniform distribution. The normal gap distribution causes 
negligibly small change in A~^. These results lead us to the conclusion that the gap anisotropy on the Fermi surface 
of MgB2 is not relevant in determining the temperature dependence of the penetration depth. 

IV. SUMMARY 

The magnetic penetration depth A(T) of a high-quality, c-axis oriented MgB2 film was obtained from mutual- 
inductance measurements in the linear-response regime. The exponential temperature dependence of A^^(r) at low 
temperatures suggests a nodeless gap on the Fermi surface. However, the data could not be described by the s-wave 
theory assuming a single gap even at low temperatures. On the other hand, the data were successfully described 
by the fuU calculation of \~'^{T) with two distinct gap values: As(0) = 2.61 ± 0.41 meV and Al(0) = 6.50 ± 0.33 
meV. At T = 0, the contribution of the small gap to the supcrfluid density was found to be 21%. Finally, two 
phenomenological models to account for gap-size distribution on the Fermi surface were considered. It was found 
that gap-size distribution in MgB2 docs not play a significant role in determining the temperature dependence of the 
penetration depth. 
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TABLE I. Summary of previously reported superconducting gap values of MgB2 superconductor. 



As (A), meV 


Al, meV 


Tool 


Refs. 


2.61 ± 0.41 


6.5 ±0.33 


Penetration depth 


this work 


1 ~ 3 


6.5 ~ 7.5 


First principle calc. 


Ref.[5] 


2.5 




Tunneling 


Ref.[8] 


3.8 


7.8 


STM 


Ref.[9] 


1.7 


7 


Point-contact spec. 


Ref.[10] 


2.8 


7 


Point-contact spec. 


Ref.[ll] 


2.45 ± 0.15 


7.0 ± 0.45 


Point-contact spec. 


Ref.[13] 


1.7 


5.6 


Photoemission spec. 


Ref.[19] 


2.7 


6.2 


Raman spec. 


Ref.[20] 


2.8 ± 0.4 




Penetration depth 


Ref.[23] 


2.61 




Penetration depth 


Rcf.[25] 
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FIG. 1. Representative complex mutual inductance M{T) of MgB2 measured using two-coil method. Inset (left) shows 
schematic diagram of measurement configuration. Inset (right) shows A^^(r) curves extracted from M{T) measured at different 
current levels. Upper cu 




FIG. 2. A"^(T) of MgB2 film at low temperatures. Thin solid and dashed lines denote best fits of A"^(0)[1 - cexp(-D/r)] 
and A"^(0)[1 - [T/Tof] to the data, respectively. 
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FIG. 3. A~^(r) of MgB2 film at low temperatures. Thin solid line denotes a best fit of Eq.(|). From this fit, A(0) = 4.29 
meV {A{0)/kBTc = 1. 31) was obtained. Inset: A^^(T) for temperatures below Tc- Thin solid line is a full BCS calculation 
assuming A(0) = 4.29 n 
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FIG. 4. A"^(T) fitted by Eq.(||). In this fit, the values of two distinct gaps As(0) = 2.57 meV and Al{0) = 5.82 meV were 
obtained. 
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FIG. 5. A^^(r) for whole temperature range below Tc. Thin solid line is a full BCS calculation of A^'^(r) assuming two 
distinct gaps. From this, refined gaps values As(0) = 2.61 meV and A_t(0) = 6.50 meV were obtained. Inset: Theoretical 
curves of A~^(r). Upper and lower curves are the contributions of the large (Al(0) = 6.50 meV) and the small (As(0) = 2.61 
meV) gaps, respectively. 
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